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1 Introduction 

The ability of an MEG system to resolve small or 
deep sources is solely determined by the achieved 
signal-to-noise ratio (SNR). Ideal MEG sensors 
should exhibit a large gain relative to the brain sig¬ 
nals and should also effectively attenuate the envi¬ 
ronmental and brain noise. Designing the system for 
large signal strength without regard for noise at¬ 
tenuation, or designing for excellent noise attenua¬ 
tion without paying attention to the signal strength 
does not produce the best results. The simultaneous 
satisfaction of large signal strength and good noise 
attenuation is not feasible in practice and the best 
sensor system represents a compromise that maxi¬ 
mizes the SNR. The paper will examine optimal and 
sub-optimal sensor designs and will illustrate their 
behavior using a synthetic signal embedded in the 
measured environmental and white noise. 

The discussion will be focused on MEG systems 
based on radial magnetometers or gradiometers, and 
noise cancellation using references, as in Fig.l.a [1]. 
In such systems, the primary sensors are located 
close to the scalp surface and are exposed to the 
brain signal and the enviro nm ental noise. The refer¬ 
ences are positioned farther from the scalp surface 
and detect mostly enviro nm ental noise. The primary 
sensors can be radial gradiometers or magnetome¬ 
ters and are both associated with a “baseline”. The 
baseline for the magnetometers enters the picture 
because even in conventional shielded rooms, the 
environmental noise is too large and the magne¬ 
tometers must be operated with some noise cancel¬ 
lation method [2], The magnetometer baseline is 
then a distance between the magnetometer primary 
sensor and the references and is usually long and 
varies widely over the scalp surface because the 
helmet dimensions are large. The primary gradi- 
ometer baselines are typically short and are identical 
for all sensors. For both types of primary sensors, 
the references are used to either synthesize higher- 
order gradiometers or adaptively subtract the noise. 
The optimum primary baseline length is determined 
by the interplay between the magnitudes of the de¬ 
tected brain signal and the noise [3], The non-bio- 
logical noise acting on an MEG system consists of 


white instrumental and low frequency environmental 
components, Fig.l.b. The important parameters are 
the white noise level, v w , the onset of low frequency 
noise, f 0 b, and the low frequency noise log-log slope, 
k, defined via Vi ow — A/f\ where A is a constant. The 
constant A is proportional to the baseline length, A 
= A 0 b/b 0 , where A and A 0 correspond to the base¬ 
lines b and b 0 , respectively. Thus if the gradiometer 
baseline is long, the magnitude of the detected envi¬ 
ronmental noise will be large and the f 0 b will also be 
large. 

The contribution of the low frequency noise depends 
on the lowest frequency of interest, fi, and the 
bandwidth of the measurement, Afi as in Fig.l.b. If 
fi is small, low frequency noise plays an important 
role, while if it is large, the low frequency noise may 
be negligible. For a constant fi, the magnitude of the 
detected environmental noise increases with in¬ 
creasing baseline b (or increasing f ob ) [3], Similarly, 
the detected brain signal magnitude also increases 
with increasing baseline [3], Since the functional 
dependencies on the baseline of the detected noise 
and the brain signal are different, the SNR (which is 
their ratio) exhibits a maximum as a function of the 
baseline. The SNR maximum defines the optimum 
operating point for the MEG system, Fig.2.a. The 
optimum baselines are typically short, in the range 
from about 2 to 8 cm and they depend on the envi¬ 
ronmental noise at a particular site [3], 

Typical noise parameters encountered within shiel- 



Figure 1: The environmental noise detected by gra¬ 
diometers depends on gradiometer baseline, (a) 
MEG system where primary sensors are radial gra¬ 
diometers or magnetometers; (b) Schematic repre¬ 
sentation of the environmental and instrumental 
noise, f Q b is the onset of the low frequency noise 
measured with gradiometers with baseline b. 
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Figure 2: (a) Illustration of the signal-to-noise ratio 
dependence on the baseline length. Sensor shell ra¬ 
dius r = 11 cm, n w = 5 jTNHz, N ave = l,k= 1, b g = 
5 cm, f = 0.1 Hz, Af= 100 Hz, dipole distance from 
the model sphere center a = 6 cm, q = 10 nA.m. 
Solid lines - SNR for various low frequency onsets, 
dashed line — optimum baseline length; (b) Exam¬ 
ples of the environmental noise parameters meas¬ 
ured within the shielded rooms at various MEG in¬ 
stallations with l s '-order radial gradiometers with 
b g = 5 cm. The blocks indicate the range of pa¬ 
rameters over all sensors of the whole cortex MEG 
system, Tu-Tuebingen, O-Osaka (night), S-Sendai, 
A-Amsterdam, To-Tokyo, P-Paris, V-Vienna, Tr-To- 
ronto. The parameters change with time, the results 
correspond to one or several measurements only. 

ded rooms at various MEG installations are summa¬ 
rized in Fig.2.b. The onset of low frequency noise 
for l st -order gradiometers range from < 1 Hz to up 
to 10 Hz, and the low frequency noise slopes k are 
in the range from about 1.5 to 3. 

To illustrate the baseline length effect on the MEG 
detection, two systems will be compared: a system 
with radial gradiometer primary sensors with near 
optimum baseline of b g = 5 cm and equipped with 
3 rd -order synthetic gradiometer noise cancellation 
[1], and a hypothetical magnetometer-based system 
with remote reference noise cancellation with aver¬ 
age baseline of b m = 26.5 cm. It will be assumed that 
both systems are subjected to the same environ¬ 
mental noise (measured in the shielded room in 
Vienna, Fig.2) and the same synthetic signals 
(spikes representative of typical MEG paradigms). 
The comparison will be performed in the limit of a 
deep dipole where the magnetometers are favored. 
In spite of that, it will be shown that the 
gradiometer-based systems provide better resolution 
even for deep sources and that the correct choice of 
the gradiometer baseline and noise cancellation are 
more important than the intrinsic signal response of 
noiseless primary sensors. 

2 Methods 

The environmental noise was measured within the 
shielded room in Vienna using l st -order radial gra¬ 


diometers with baseline b g = 5 cm and synthetic 3 rd - 
order gradiometers-based on them [4], The noise for 
both sensor types was filtered to the required fre¬ 
quency range of fi to fi + Af using 10 th -order But- 
terworth filters. The rms noise value, v rms , in the 
bandwidth Af was determined as a function of fi for 
fi in the range from 0.1 to 20 Hz. 

The measured rms noise, Vrm S (l st -order gradiometer, 
b g = 5 cm) was then scaled by the baseline ratio to 
estimate the noise of a magnetometer with b m = 26.5 
cm. At low fis, the Vrm S is dominated by the low fre¬ 
quency noise and the scaling by the baseline ratio is 
straightforward. Similarly, at high fis, the v rms is 
dominated by white noise (the low frequency com¬ 
ponent is absent) and the noise is independent of the 
baseline. At intermediate fis, the noise is a mixture 
of white and low frequency components. In this 
range, the magnetometer rms noise was estimated by 
extrapolation of both components and adding them 
together in an rms fashion. This extrapolation was 
used only for the SNR estimation; the time plots 
presented in Fig.4 were obtained for fi either in the 
low frequency or white noise regimes. 

The synthetic signal was constructed as a Gaussian 
spike with variable width. For the magnetometer- 
based system, the amplitude of the spike was se¬ 
lected to be 3 pT and the spike was added to the un¬ 
averaged noise (N av e = !)• This signal is characteris¬ 
tic of an interictal spike. Alternatively, the signal 
amplitude was set to 300 fT and the noise was re¬ 
duced by averaging (N ave = 100), to mimic an 
evoked response. In both cases, however, identical 
SNRs are obtained. The magnitude of the detected 
brain signal depends on the source depth and the de¬ 
vices with shorter baselines will detect smaller sig¬ 
nals. The ratio of gradiometer-based (b g = 5 cm) to 
magnetometer-based (b m = 26.5 cm) brain signal 
amplitudes was set to 0.7, corresponding to the lim¬ 
iting behavior when the source approaches the 
model sphere center (note that in this limit, the de¬ 
tected signal magnitudes will approach zero for all 
devices, however, the magnetometer-to-gradiometer 
amplitude ratio is 0.7). This choice of the signal am¬ 
plitudes provides the greatest possible advantage for 
the magnetometer-based systems. 

The signal width was varied from 20 msec (repre¬ 
senting early evoked responses) to more than 1 sec 
(representing voluntary finger motion, CNV, etc.) 
and the signal was subjected to the same filtering as 
the noise. 

The resulting noise and signal were used to either 
estimate the dipole resolution Aq, or directly calcu¬ 
late the SNR. The dipole resolution was determined 
from a 99% probability of finding the dipole within 















a 1 cm edged cube [5], The SNR = s pp /(2V2v rms ), 
where s pp denotes the peak-to-peak signal and Vrm S is 
the rms value of the noise, both after filtering. Both 
the dipole resolution and the SNR were studied as a 
function of the minimum frequency of interest, ft. 

3 Results 

The dipole resolution and the SNR as a function of 
the minimum frequency of interest are shown in 
Fig.3.a for a synthetic 3 rd -order gradiometer, a mag¬ 
netometer-based system, and a pure magnetometer. 
When ft < 7 Hz, the best dipole resolution is ob¬ 
tained for the 3 rd -order gradiometer (because of its 
efficient noise cancellation). Only at higher frequen¬ 
cies do the magnetometer-based systems become 
slightly better, because they detect larger signals (a 
factor of 1/0.7 for a dipole at the sphere center). The 
pure magnetometer resolution is the worst in the 
whole investigated range of fts. 



Figure 3: Comparison of synthetic 3 rd -order gradi¬ 
ometer- and magnetometer-based system perform¬ 
ance as a function of the minimum frequency of in¬ 
terest, fa. The noise was measured within a shielded 
room in Vienna with a 143 channel MEG system 
equipped with the synthetic 3 rd -order gradiometer 
noise cancellation [4], The noise corresponds to the 
channel MLC34, a white noise level v w = 5 fT/fHz, 
ratio of the gradiometer to the magnetometer signal 
strengths = 0.7 (the limit of a dipole at the model 
sphere center), primary gradiometer baseline b g = 5 
cm, magnetometer-based system baseline b m = 26.5 
cm (ratio b m /b g = 5.3). Dashed thick lines - magne¬ 
tometer-based system, thin dashed line — pure mag¬ 
netometer, solid lines - synthetic 3 rd -order gradi¬ 
ometer. (a) Equivalent current dipole resolution 
[5J; (b) a/SNR for signal width of 50 msec. 


The SNR for one channel is shown in Fig.3.b. 
Again, for ft < 7 Hz, the 3 rd -order gradiometer SNR 
is better than that of the magnetometer-based system 
and for ft > 7 Hz, the magnetometers become 
slightly better, because they detect larger signals. 

The time traces for the magnetometer- and gradi- 
ometer-based systems with ft = 1 and 10 Hz are 
shown in Fig.4. For ft = 1 Hz (well within the low 
frequency noise regime), the magnetometer barely 
resolves the signal and its SNR is a marginal 1.8, 
while the 3 rd -order gradiometer exhibits a large SNR 
of about 56. For ft = 10 Hz (in the white noise re¬ 
gime), both instruments resolve the signal clearly, 
however the magnetometer SNR of 27.5 is larger 
than the gradiometer SNR of 19.1. At high ft, most 
of the noise has been filtered out and the difference 
is due to the magnetometer detecting a larger signal. 
If the magnetometer to gradiometer signal magni¬ 
tude ratio is based on a more superficial dipole, the 
difference is correspondingly smaller (e.g. for a = 6 
cm the gradiometer SNR is 23.2 instead of 19.1). 

The SNR depends on the signal frequency content, 
i.e. the signal width. This is illustrated in Fig.5 for 
signal widths of 20 and 100 msec. In both cases, for 
ft < 7 Hz, the 3 rd -order gradiometer SNR is larger 



Figure 4: Comparison of the time traces for a syn¬ 
thetic 3 rd -order gradiometer and a radial magne¬ 
tometer for the minimum frequencies of interest of fa 
= 1 and 10 Hz. The environmental noise was meas¬ 
ured within the shielded room in Vienna (as in 
Fig.2.b) and was filtered with 10 lh -order low and 
high pass Butterworth filters to the bandwidths of 1 
to 21 Hz and 10 to 30 Hz. The signal was a Gaus¬ 
sian spike 50 msec wide with 3 pT amplitude and 
Nave = 1 (or 300 fT amplitude and N ave = 100). The 
SNR is indicated by the number associated with 
each plot, (a, b) The filtered signal for f = 1 Hz; (c, 
d) The filtered signal for fa = 10 Hz; (a, c) Magne¬ 
tometer-based system; (b, d) Gradiometer-based 
synthetic 3 rd -order gradiometer. 























Figure 5: Square root of the one-channel SNR in the 
presence of measured noise as a function of the 
minimum frequency of interest, j), for magnetome¬ 
ter- and gradiometer-based systems. Noise and sig¬ 
nal amplitudes are as described in captions for 
Figs.3 and 4. Solid line - gradiometer-based system 
with synthetic 3 rd -order gradiometer noise cancel¬ 
lation, dashed line — magnetometer-based system. 
Ratio of the gradiometer to the magnetometer signal 
strengths = 0.7 (limit of the dipole approaching the 
model sphere center), (a) Signal width 20 msec, 
characteristic of early evoked responses; (b) Signal 
width 100 msec, characteristic of late evoked re¬ 
sponses. 


than that of the magnetometer. For the 20 msec 
width in Fig.5.a, the magnetometer SNR above fi = 
7 Hz increases roughly to that obtained with the 3 rd - 
order gradiometers at lower frequencies. For the 100 
msec width in Fig.5.b, the magnetometer SNR re¬ 
mains low even at high frequencies and it never at¬ 
tains the large magnitudes observed for the 3 rd -order 
gradiometer SNR. 

4 Discussion 

It was shown that the best source resolution, or the 
ability to detect deep sources, depends on the SNR. 
Two systems were investigated: First, a magne¬ 
tometer-based system with noise cancellation and a 
long effective baseline (26.5 cm). This system ex¬ 
hibits a large gain relative to the brain sources, how¬ 
ever, its baseline is too long and the system has ex¬ 
cessive sensitivity to noise. Second, a l st -order gra¬ 
diometer-based system with a short optimized base¬ 
line (5 cm) and with synthetic 3 rd -order gradiometer 
noise cancellation. This system has a smaller gain 
relative to the brain sources (especially when they 
are deep); however, it exhibits efficient noise can¬ 
cellation. The results indicate that even for deep 
sources, where the magnetometer-based system de¬ 
tects larger signals, optimized gradiometers have 
larger SNR and better ability to resolve smaller 


The characteristics of the two systems were exam¬ 
ined by adding synthetic signals to the measured en¬ 
vironmental noise and then evaluating the SNR and 
the magnitude of a minimum resolvable dipole. The 
performance is a function of the signal frequency 
content and the minimum frequency of interest, fi. 
For lower fis, the SNR and the ability to resolve di¬ 
poles are vastly better for the optimized gradiometer 
system than for the magnetometer system. For larger 
fis, the magnetometers become slightly better, how¬ 
ever, even then their SNR is typically smaller than 
the gradiometer SNR for smaller fis. Thus, if each 
device was operated at its optimum fi, the optimized 
gradiometer performance would be much better than 
that of the magnetometers and the gradiometers 
would resolve smaller sources at any depth. If the 
gradiometers were forced to operate at the optimum 
fi for the magnetometers, then the magnetometer 
SNR would be slightly better. However, if the mag¬ 
netometers were forced to operate at the optimum fi 
for the gradiometers, their SNR would be exceed¬ 
ingly low and no useful signal would be recovered. 
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